Shear-induced Diffusion in Dilute Suspensions of Charged Colloids 
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We propose a model for the nonequilibrium enhancement of colloidal self-diffusion in an externally 
imposed shear flow in charged systems. The enhancement of diffusion is calculated in terms of the 
electrostatic, two-body interactions between the particles in the shear flow. In the high shear 
rate, low volume fraction limit in which our model is valid, we compare these calculations to the 
experiments of Qiu et al. [Phys. Rev. Lett. 61, 2554 (1988)]. 
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Colloidal suspensions ||, |) serve as a unique model sys- 
tem for both atomic liquids and solids. Due to their in- 
herently longer length and time scales, colloidal systems 
offer the opportunity to probe the many-body, nonequi- 
librium dynamics of interacting systems through the ap- 
plication of e.g. moderate shear. Such issues involving 
systems driven far from equilibrium remain among the 
principal puzzles of modern statistical mechanics. 

In this letter we discuss the origin of a particularly 
intriguing and apparently generic feature of particulate 
suspensions driven out of equilibrium by the application 
of shear flow. Since its initial observation by Qiu et al. 
H in a charged colloidal suspension, it is now well-known 
that the effective single-particle diffusion constant grows 
under shear. Similar enhancement of self-diffusion under 
applied shear been observed in concentrated non-colloidal 
(i.e. non-Brownian or high Peclet number) suspensions 
both experimentally [|[ [| |) and in simulations Q. It 
should be pointed out that the diffusion enhancement 
under consideration occurs in the gradient and vorticity 
directions, and is thus unrelated to the better understood 
Taylor dispersion which contributes to the effective dif- 
fusion along the flow direction. 

Since laminar, zero Reynolds number flows obey time 
reversal symmetry, one is forced to look elsewhere for the 
symmetry breaking required for understanding the phe- 
nomenon of nonequilibrium diffusion enhancement. In 
uncharged, non-colloidal suspensions it has been shown 
numerically that three-particle hydrodynamic interac- 
tions are chaotic || [|. These chaotic interactions have 
been previously implicated as a source of non-thermal 
noise in nonequilibrium suspensions p0| . In addition, fi- 
nite particle roughness |Tl| can break the symmetry of 
low Reynolds number hydrodynamics effectively by cut- 
ting off lubrication forces at small interparticle separa- 
tions. 

In interacting colloidal systems there is, however, an- 
other source of symmetry breaking due to the elec- 
trostatic interactions between the colloids [l^] (See 
Fig. [I]) . The importance of the electrostatic interactions 
in the diffusion enhancement observed by Qiu et al. is 
demonstrated by the fact that a decrease in the Debye 



screening length suppresses the effect. We focus exclu- 
sively on this experiment and the role of the electrostatic 
interactions since these interactions are more analytically 
tractable than the chaotic, many-body hydrodynamics 
that presumably generates a similar diffusion enhance- 
ment in non-Brownian hard-sphere suspensions. 

Our calculation proceeds as follows: We consider the 
randomization of the trajectory of a particle, which we 
refer to as the scattering center, due to its electrostatic 
interaction with a spatially random collection of identical 
particles that are carried past it in the shear flow. This 
is done in two steps. We first compute the trajectory of 
a charged colloid as the macroscopic shear flow carries 
it past the scattering center. From this calculation we 
determine the displacement of the scattering center itself. 

Given a spatially random ensemble of such particles 
being carried by the scattering center by the simple shear 
flow (uniform rate of shear equal to 7), we compute the 
second moment of the displacements of the scattering 
center due to interactions with the particles flowing past 
it. The product of this second moment with the rate 
of scattering events (~ 7) gives the shear enhancement 
of the diffusion constant. The result of this calculation 
is the diffusion enhancement: AD ~ |7|". We compare 
both the magnitude of this diffusion enhancement and its 
dependence upon shear rate (i.e. a) with experiment. 

Qualitatively, the trajectory of a charged particle un- 
der the combined influence of the macroscopic shear flow 
and the electrostatic interaction with the scattering cen- 
ter can be described as follows: The shear flow brings a 
particle toward the scattering center along a given stream 
line. While the particle is within a few screening lengths 
of the scattering center, the electrostatic interaction dis- 
places the particle from its stream line. However, once 
the particle has been carried by the combination of flow 
and electrostatic interaction to a distance of more than a 
few screening lengths, the particle resumes its trajectory 
along a different streamline than the one it entered on. 
See Fig. ^ for computed trajectories. 

There are two important physical limits which result 
in qualitatively different particle trajectories and a dif- 
ferent dependence of diffusion enhancement upon shear 
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rate. The limits correspond to the dominance of one of 
two independent forces controlling the particle trajecto- 
ries, hydrodynamic drag and the interparticle, electro- 
static interaction. 

In the case where the electrostatic interaction is arbi- 
trarily strong so that it completely controls the motion 
of a particle in the interaction zone, we recover a sim- 
ple linear dependence of the diffusion enhancement with 
shear. A particle enters the interaction zone along some 
stream line. Now the dominant electrostatic repulsion 
between that particle and the scattering center drives 
the particle out of the interaction zone. It then resumes 
its straight-line motion along a different streamline. The 
displacement of the particle is thus on the order of the 
interaction zone radius and independent of the imposed 
shear rate. The scattering rate, however, is proportional 
to the shear rate, so one the diffusion enhancement scales 
linearly with shear rate, AD ~ |7| in the low shear rate 
regime where the particle-particle interaction is entirely 
dominated by the electrostatics. See Fig. 2a. 

On the other hand, if the electrostatic interaction is 
weak enough, or, equivalently, if the shear rate is high 
enough, a particle will be carried predominantly by the 
shear flow through the electrostatic interaction zone sur- 
rounding the scattering center. In this case the residence 
time of the particle in the interaction zone is inversely 
proportional to the shear rate. In the low Reynolds num- 
ber limit (which we always assume) the displacement of 
the particle from its initial stream line is proportional to 
the time integral of the force acting on it. In this case 
that time integral will be proportional to 7~ 1 . Since the 
rate of scattering events will still be proportional to 7, 
we expect that the enhancement of diffusion will plateau 
at high shear rates. See Fig. 2b. 

The above, heuristic cases describe only the limiting 
cases of strong and weak shear. Not surprisingly, for 
physical values of the ratio of these two forces including 
those encountered in experiment, the predicted result for 
the shear enhancement interpolates between these two 
limits. We find for a broad range of shear rates and 
for physically relevant electrostatic interaction parame- 
ters that the shear enhancement of diffusion scales as 
AD~\j\ - 7 . 

Our calculation cannot be extended to zero 7. We esti- 
mate a lower shear-rate cutoff for our analysis by noting 
that we assume that particles follow stream lines towards 
and away from the scattering center except when they are 
within the interaction zone surrounding each particle. If 
the shear rate is so low that particles can diffuse through 
the interaction zone, the deterministic trajectory calcula- 
tions are no longer valid. The cross-over from our model 
calculations to quadratic scaling occurs at a critical shear 
rate 7* ~ k 2 Do where Dq is the Brownian diffusion con- 
stant. At this cross-over shear rate a particle diffuses 
a Debye length in the same time as it is advected that 
distance by the flow. For the experimental system in 



question 7* = 10 2 s _1 . In addition, we note that the re- 
duction of many-body interactions into a series of simple 
two-body interactions can only be justified in the limit 
of small particle number density. 

This result for the scaling of the shear enhance- 
ment of diffusion with shear rate is in good agree- 
ment with noncquilibrium, Brownian-dynamics simula- 
tions jD| The experimental data ||, on the other 
hand, is consistent with a power law of unity. To resolve 
this discrepancy, we point out that our calculation is not 
applicable to the lower shear rate experimental data. In 
this lower shear rate regime, we expect that the shear- 
rate dependence of AD is stronger than linear fl6|| so that 
the effect of the cross-over of the exponent from greater 
than one to 0.7 leads to a larger apparent exponent in the 
data]!?]]. Secondly, we note that, over the single decade 
of shear-rate data available, it is problematic to distin- 
guish between our proposed exponent of 0.7 and the lin- 
ear dependence on shear rate concluded by Qiu et al. 
In addition, we will show, that by making reasonable as- 
sumptions about the electrostatic interactions present in 
the experimental system of Qiu et al, our calculation 
of the magnitude of the shear enhancement effect agrees 
well with the experiment. 

In an earlier mode-coupling study Indrani it et al. G3| 
attribute the diffusion enhancement to the suppression of 
the effective friction, (, experienced by a concentration 
fluctuation. Noting that D ~ T/£ they calculate the 
contribution to £ from the interaction of a particle con- 
centration fluctuation with all other thermally generated 
concentration fluctuations in the medium. The principal 
effect of the shear flow is to destroy these concentration 
fluctuations and thereby decrease the effective friction. 
We, on the other hand, propose that the principal effect 
of the shear flow when considered along with the electro- 
static interparticle interactions is to increase the effective 
temperature in the nonequilibrium system. 

The particle trajectories are described by the following 
differential equation fig]: 



dY 



Y- 



H Y + l 



sin (p 



(1) 



where Y — kR is the distance of the incoming particle 
from the scattering center measured in Debye lengths k _1 
and £ = cot is the polar angle. See Fig. |l| Eq. [I] ac- 
counts for the interplay of the hydrodynamic drag force 
by the fluid represented by the first term on the LHS of 
and electrostatic repulsion between particles represented 
by second term on the LHS of Eq. |l|. The relative im- 
portance of these two forces is measured by the dimcn- 
sionless parameter H, which represents the ratio between 
the electrostatic interparticle interaction and the hydro- 
dynamic drag. H is defined by 
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where e r eo is the dielectric constant of the fluid, $ a the 
apparent surface potential of the particles, and C the hy- 
drodynamic friction factor. All the details of the elec- 
trostatic interaction for a given colloidal system can be 
subsumed into the parameter C defined above. 

Eq. H neglects hydrodynamic interactions between par- 
ticles which is acceptable if particles do not make close 
passes. In a dilute suspension of highly charged colloids 
the number of such close passes should be negligible and 
thus not greatly contribute to the effective diffusion con- 
stant. This view is supported by the original experiments 
which showed that the shear-enhancement of diffusion 
disappears with sufficient screening of the electrostatic 
interaction ||. 

To determine the trajectory of a scattering particle we 
integrate Eq. |l| from initial conditions such that at large 
distances from the scattering center (i.e. x — ► — oo) 
the incoming particle has an impact parameter ro and 
has a polar angle ipo measured from the z-axis in the yz 
plane. As a result of the azimuthal symmetry about the 
x-axis, every trajectory is confined to a plane of constant 
ip. However, because y is the velocity gradient direction 
of the imposed shear flow, the rate of incoming particles 
to the scattering center will depend on ip. We show in 
Fig. U two numerically integrated trajectories. The first 
shows the result obtained in the electrostatically domi- 
nated regime (high H) and the second shows the shear 
flow dominated regime (low H). 




where P(s)ds is the probability per unit time of a collision 
to occur with impact parameter s = (fo,tpo), where fo = 
ktq. The terms in the parentheses are the displacements 
of a given particle in the respective directions. Assuming 
a random distribution of particles, P{s) is given by the 
particle flux through the area tq dro dpo ■ 

P(s)d 2 s=- — -, — ^rj sintpo rl dr d(p (5) 
47T (na) 

The integration domain in Eqs. |3j^| only accounts for par- 
ticles flowing in from the left (y > 0). The values for 
the elements of the diffusion tensor have been doubled 
to account for the contribution of the particle flux com- 
ing from the right (y < 0). As expected, the shear rate 
7 affects shear-enhanced diffusion in two distinct ways: 
through the particle flux (~ 7) and the step size, which 
depends on H = C/7 in a complicated way that we ex- 
plore below. 
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FIG. 1: Coordinate system for two-particle collision in a shear 
flow; the right particle is located at the origin of the frame of 
reference. 

We determine the total displacement of the in-coming 
particle normal to the streamlines by finding the change 
in the radial distance Af = k(t\ x >oc — ro) of the par- 
ticle's trajectory from the i-axis comparing the asymp- 
totic out-state where x — > 00 to the initial impact pa- 
rameter. By symmetry, each particle (the incoming par- 
ticle and the scattering center) moves half of the distance 
Ar in the laboratory frame. 

The effective enhancement of the diffusion constant in 
both the velocity gradient (y) and vorticity (z) directions 
can be calculated from the second moment of the particle 
displacements in the appropriate directions: 
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FIG. 2: Trajectories in the x — y— plane (i.e. (pa — tt/2) for 
the collision between two electrostatically interacting colloids 
in a shear flow; H = 10 6 in the upper figure and H — 1 in 
the lower figure; trajectories are plotted for different impact 
parameters fo. The dashed semicircles denote the volume of 
the interaction zone surrounding the scattering center. 

To explore the shear rate dependence of the diffusion 
enhancement, it is convenient to express the diffusion 
constants in the following form using Eqs. §@, and § 
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AD nn = 



(0) 



where I a is a dimensionless integral and a = y, z. The 
entire shear rate dependence of the two diffusivities is 
contained in 1/H ■ I a (H) where H ~ 1/7. For small 
values of H (large 7 or low interaction), the H de- 
pendence of I a (H) is linear as can be seen by simple 
perturbation theory in H. AD aa at these large shear 
rates becomes independent of H. For large H , however, 
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the dependence on H weakens and the curves are de- 
scribed well by I a (H) ~ H 3 so that the diffusion en- 
hancement, which scales as 1/H ■ I a (H) takes the form: 
AD a ~ 1/H ■ I a (H) ~ H-° J ~ 7 7 . That this power 
law holds over a wide range of shear rates is one of the 
central results of this letter. Such power law scaling is 
difficult to distinguish from an exponent of unity over 
the range of shear rates where we expect the model to 
be valid (100 — 1000 s _1 ) The scaling is not inconsistent 
with the experimental data of Qiu et al. In addition, the 
computed exponent agrees well with previous simulations 
of the system Q. 

The dependence of the diffusivities upon 1/H demon- 
strating the cross-over from the 0.7 power law to the H 
independence is shown in Fig. ||. This scaling of the dif- 
fusivity with shear rate is distinct from that observed in 
non-colloidal suspensions, where shear-induced diffusion 
grows monotonically with j. At small shear rates the dif- 
fusion enhancement is anisotropic: AD yy / AD ZZ w 1.7. 
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FIG. 3: Scaling of diffusion with shear rate represented by 
the equivalent plot of I(H)/H vs. 1/H. 

The calculated magnitude AD is also consistent with 
the experimental data. To test this point, we need to find 
C . It is determined from a combination of the electro- 
static properties of the suspension including the surface 
potential of the colloids $ a , the Debye screening length 
and the particle radius as well as the inverse of the 
effective particle mobility. Since we determine sphere mo- 
bility empirically from single-particle diffusion constant 
of the uncharged system, we have only to know the elec- 
trostatic parameters of the charged system in order to 
determine C, and thus determine the calculated diffusiv- 
ity enhancement for the sheared suspension. 

Unfortunately, Qiu et al. do not report a value for 
the surface potential of the charged colloids. In spite of 
the fact that the experimental system is not well enough 
characterized to fix C, we can, by using the surface po- 
tential as a fitting parameter, force H to be in the range 
where we find the scaling relation AD ~ |7| ' 7 . To 
achieve this we set $ a = 200mV, which is not atypi- 
cal in such charged colloidal systems. When H is set 



to 10 5 and 7 is chosen equal to be 500 s _1 — in the 
middle of the experimental range — our model predicts 
AD ZZ ~ 3 • 10~ 12 m 2 /s. This prediction is reasonably 
close to the experimental result of AD ZZ ~ 1-10~ 12 m 2 /s, 
considering the strong dependence of AD on the electro- 
static parameters, which are known imprecisely at best 



(e.g. AD 



.-4.1 „0.6kq 



). The agreement between ex- 



periment and theory demonstrates only the plausibility 
of the latter. Diffusion measurements under shear need 
to be performed on better characterized systems to rigor- 
ously test the theory. Furthermore, systematic variation 
of the electrostatic parameters will provide insight into 
the validity of our scaling predictions. 
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